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Abstract. We employ a modified composite day extension 
of the Hocking (2005) analysis method to study gravity 
wave (GW) activity in the mesosphere and lower thermo- 
sphere using 4 meteor radars spanning latitudes from 7° S to 
53.6° S. Diurnal and semidiurnal modulations were observed 
in GW variances over all sites. Semidiurnal modulation with 
downward phase propagation was observed at lower lati- 
tudes mainly near the equinoxes. Diurnal modulations occur 
mainly near solstice and, except for the zonal component at 
Cariri (7° S), do not exhibit downward phase propagation. At 
a higher latitude (SAAMER, 53.6° S) these modulations are 
only observed in the meridional component where we can 
observe diurnal variation from March to May, and semidi- 
urnal, during January, February, October (above 88 km) and 
November. Some of these modulations with downward phase 
progression correlate well with wind shear. When the wind 
shear is well correlated with the maximum of the variances 
the diurnal tide has its largest amplitudes, i.e., near equinox. 
Correlations exhibiting variations with tidal phases suggest 
significant GW-tidal interactions that have different charac- 
ters depending on the tidal components and possible mean 
wind shears. Modulations that do not exhibit phase variations 
could be indicative of diurnal variations in GW sources. 

Keywords. Meteorology and atmospheric dynamics (Waves 
and tides) 


1 Introduction 

Large-scale dynamics in the mesosphere and lower ther- 

mosphere (MLT) are strongly influenced by smaller-scale 
gravity waves (GWs) owing to significant GW momentum 


transport and deposition accompanying GW propagation and 
dissipation. Because GW dissipation is strongly influenced 
by large-scale wind shears, GW propagation and momen- 
tum transport are expected to be strongly modulated by 
the various tidal and planetary wave (PW) motions (e.g., 
Walterscheid, 1981; Fritts and Alexander, 2003). Indeed, a 
number of previous observations have revealed strong mod- 
ulations of GW variances and momentum fluxes, suggest- 
ing significant interactions with the tides and PWs. MF and 
VHF radars have indicated (1) peaks in GW variance spec- 
tra at tidal and PW periods (Isler and Fritts, 1996; Manson 
et ah, 1998) and (2) an approximate anti-phase relation- 
ship between tidal winds and GW momentum fluxes (Fritts 
and Vincent, 1987; Wang and Fritts, 1991). Modelling stud- 
ies have suggested both tidal amplitude reductions and/or 
phase advances (Forbes et ah, 1991; Lu and Fritts, 1993; 
Meyer, 1999) and amplitude enhancements (Mayr et ah, 
1998), depending on the GW parameterization used. Ortland 
and Alexander (2006) demonstrated that these influences de- 
pend strongly on the phase speed distribution and isotropy of 
the GW spectrum. 

Meteor radars can also contribute to studies of GW-tidal 
and GW-PW interactions through definition of the large- 
scale motions and the GW variances and momentum fluxes. 
Studies of GW variances or momentum fluxes to date have 
addressed primarily seasonal variations (e.g., Antonita et ah, 
2008; Beldon and Mitchell, 2010; Fritts et ah, 2010, 2012a; 
Placke et ah, 2011a), most using the Hocking (2005) anal- 
ysis method or a derivative thereof. The Hocking (2005) 
method in principle allows the use of meteor radar data to 
estimate the large-scale motions as well as all GW variances 
and momentum fluxes. This has enabled significant advances 
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Fig. 1. Available data series for the meteor radar located at Cariri. The blank blocks indicate data gaps. 
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Fig. 2. The same as Fig. 1 but for the meteor radar at CP. 


in studies of GW propagation and effects in the MLT, given 
the worldwide network of meteor radars that operate almost 
continuously. 

Despite the benefits of the Hocking (2005) method, there 
are biases that arise in the estimates of GW variances and 
momentum fluxes due to the spatial and temporal averag- 
ing required to form these estimates with sufficiently small 
uncertainties. Vertical wind shears and temporal wind vari- 
ations both contribute to increased apparent GW variances 
and biased momentum fluxes unless corrections are applied. 
Placke et al. (2011a, b) attempted to account for vertical wind 
shears, while Andrioli et al. (2013) examined both effects and 
provided a modified composite day (MCD) analysis method 
that overcomes the more general problem. 

In this paper we present GW variance estimates obtained 
by applying the MCD version of the Hocking (2005) tech- 
nique using meteor radars at several latitudes. Findings in- 
clude diurnal and semidiurnal modulations of GW variances 
exhibiting good correlations with tidal wind shears in several 
cases. Our findings provide additional evidence for strong 
GW-tidal interactions observed in previous studies. These 
dynamics probably also include generation of additional 


GWs accompanying these interactions, given the expectation 
for secondary GW generation in regions of strong dissipation 
and momentum deposition (Vadas and Fritts, 2001, 2002). 
Meteor radar measurements are probably not able to mea- 
sure these responses, however, because any secondary GWs 
generated at these altitudes will only yield large amplitudes 
and variances at higher altitudes. 

2 Method 

Hocking’s (2005) analysis assumes that winds in the MLT 
region are uniform on a horizontal plane and that any dif- 
ference between a given measured meteor radial velocity 
and the fitted radial wind velocity is due to the contribu- 
tion of GWs. On this basis it is possible to compute the 
meridional, zonal, and vertical fluctuating wind velocities 
and the vertical flux of horizontal momentum (see Hocking, 
2005, for details). In the present paper we apply the MCD 
extension of the Hocking (2005) method to data from four 
all-sky meteor radars. The MCD analysis involves a simple 
change in the way in which traditional composite day esti- 
mates are constructed. This involves a pre-analysis that infers 
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Fig. 3. The same as Fig. 1 but for the meteor radar at SM. 



Fig. 4. Daily meteor distributions, for 1 June 2005, illustrating the 
azimuth patterns of the three Brazilian radars employed in this 
work. (Left) Meteor distribution including all the zenith angles and 
(right) only the meteors used in the analysis. The location of each 
of the radars is listed in each panel. 


the differences between individual radial velocities and best 
fit radial velocities (hereafter referred to as u' ad ) for each time 
interval and day separately, and associates them with the me- 
teor position information throughout the entire data set. After 
this procedure, we then apply the Hocking (2005) technique 
and compute the GW wind variances and vertical fluxes of 
horizontal momentum. In this way we can use a composite 
day analysis but avoid the effects of day-to-day changes in 
tides and PWs. Furthermore, using MCD we can accumulate 
more echoes in each height/time interval, and the larger the 
meteor count the greater is the confidence in the results. 

The data were analysed using a height interval of 4 km 
width (centered at 82, 85, 88, 91, 94, and 98 km), a 3 h time 
interval (centered at 01:00, 03:00, 05:00, 07:00, 09:00, 11:00, 


13:00, 15:00, 17:00, 19:00,21:00, 23:00 UT), and zenith an- 
gles between 15 and 50 degrees. The latter constraint avoids 
spurious contributions to large GW variances due to large 
vertical velocities at small zenith angles and range errors 
due to zenith angle uncertainties at low elevation angles. 
This analysis allows us to study GWs with periods less than 
~3h, vertical wavelengths less than ~ 5-10 km, and hori- 
zontal wavelength less than ~ 1 80 km. Additional details of 
the method can be found in Andrioli et al. (2013). 

We analyze the data from three SKiYMET meteor 
radars and the Southern Argentina Agile MEteor Radar 
(SAAMER), which are well distributed in latitude, including 
Sao Joao do Cariri (Cariri) at 7° S; 36° W, Cachoeira Paulista 
(CP) at 23° S; 45° W, Santa Maria (SM) at 30° S; 54° W, and 
SAAMER on Tierra del Fuego at 53.8° S; 67.8° W. With the 
exception of SAAMER, where we used only the data for 
2010, our analysis employs all available data. The results 
presented here represent monthly means for each radar for all 
available years. Figures 1, 2, and 3 show the available data se- 
ries from Cariri, CP, and SM, respectively. The blank blocks 
represent data gaps when, for various reasons, the radars did 
not operate. The SAAMER data used are from January to De- 
cember of 2010 when the radar worked nearly continuously, 
with only five days without data. We can see from these fig- 
ures that there are approximately four, nine, and five years of 
data for Cariri, CP, and SM, respectively. 

As the Hocking (2005) analysis is a generalization of the 
Vincent and Reid (1983) dual-beam technique, the azimuthal 
meteor distribution has an important influence on the accu- 
racy of the analysis. Fritts et al. (2012b) have shown that 
the best efficiency of Hooking’s analysis is achieved when 
the meteor counts in opposing directions are large. Figure 4 
shows the daily meteor distributions in zenith angle and az- 
imuth for the three Brazilian radars observed on 1 June 2005. 
Note that SM uses a crossed Yagi transmitting antenna and, 
for this reason, it has a more uniform radiation pattern. In 
contrast, CP and Cariri use linear Yagi transmitting anten- 
nas. Nevertheless, for the zenith angles from 15° to 50° em- 
ployed for our studies, all radars have reasonably uniform 
azimuthal sensitivities. This means that the MCD version of 
the Hocking (2005) method should work well for all radars. 
Daily meteor distributions for SAAMER are shown in Fritts 
et al. (2010, 2012a, b). 
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Fig. 5. Zonal component of variances (left side) and the total wind (right side) averaged from 2005 to 2008 over Sao Joao do Cariri. All 
panels show the height/time variation of their respective variable. The altitude is measured in km, and the x axis represents Universal Time 
(UT) beginning at midnight, and each line of the variances contour plot corresponds to 50 m 2 s~ 2 . The total wind is measured in m s~ 1 and 
the white dashed line indicates where the values are zero. 


It should also be remembered that GW variances are likely 
to be somewhat more uncertain where meteor counts are 
lower, hence at higher and lowest altitudes. This is because 
Hocking’s analysis is not as well conditioned with fewer me- 
teor counts. 

2.1 Method for removing tidal biases from the Hocking 
variance analysis 

Andrioli et al. (2013) developed an empirical technique to 
remove the apparent variances from the Hocking analysis, 
leaving only estimates of the radial velocities due to GW mo- 
tions (e.g., the MCD method). In the present work we use the 
MCD method in order to analyse the GW variances and their 
correlations with tidal winds and shears over the 4 radars. 
This method involves three steps. Step 1 infers the tidal fields 
and total variances using Hocking’s method. Step 2 employs 


the tidal parameters obtained in Step 1 as input for a fit to 
the large-scale wind field that varies smoothly in space and 
time. This allows, in Step 3, more accurate estimates of indi- 
vidual radial velocities and resulting GW variances. The fits 
in Step 2 are performed using the measured meteor parame- 
ters, including the position of each detected meteor, and re- 
placing the measured radial velocity by that calculated from 
the model, using Eqs. (1), (2), (3) and (4) below. Finally, we 
subtract the apparent variances from the total, leaving only 
the variances due to GWs. Additional details are provided by 
Andrioli et al. (2013). 

U(x,y, z, t) = U M + Ub(z, t) sin (2n(t - <5 UD )/7b) 

+Usd(z, f)sin(2^(t -<5usd)/7sd) (1) 

V(x,y,z,t ) = Vm- Vd(z, Ocos (2n(t — <5vd)/7d) 

- Vsd(z, f)cos(27r(f - <5 V sd/)/7sd) (2) 
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Fig. 6. The same as Fig. 5, but for the meridional component. 


3 Results 


W(x, y, z, t) = 0 (3) 


V m d = U(x,y,z, 0 sin# cos + V{x,y,z, Osin# sin <p 

+ W(x,y,z,t)cos6 (4) 

Here, U, V, and W , are zonal, meridional and vertical wind 
velocities, V r ad is the radial velocity calculated for each me- 
teor zenith and azimuth position, #, (p , respectively, sub- 
scripts M, D, and, SD denote mean wind and diurnal and 
semidiurnal tides (with amplitudes assumed to vary smoothly 
in space and time), {x,y,z) is the meteor position, and t 
is the time when each meteor was detected. (<5ud, <5vd ) and 
(<5usd, <5vsd) are diurnal and semidiurnal tidal phases vary- 
ing according to the inferred vertical wavelengths, subscripts 
U and V refers respectively to zonal and meridional compo- 
nent, and 7’d and'/sn are the respective tidal periods, 24 and 
12 h. 


Several studies have been performed analyzing data from the 
three radars at Cariri, CP, and SM with a focus on tides and 
PWs (Andrioli et ah, 2009; Batista et ah, 2004; Lima et ah, 
2004, 2005, 2006, 2007; and Buriti et ah, 2008). Although 
Clemesha and Batista (2008) and Clemesha et ah (2009) pre- 
sented some studies of wind variances related to GWs using 
these data, the possible effects of tidal contamination were 
not taken into account. As noted above, Andrioli et ah (2013) 
showed that GW variances are typically contaminated by 
tides and suggested a method for removing this contamina- 
tion, thus allowing their accurate estimate. For this reason, 
we re-analyze all previous data using our new MCD method. 
We also analyze one year of SAAMER data in order to ex- 
tend our analysis to higher southern latitudes. 

Figure 5 shows the zonal variances (left side) and total 
wind (right side) for each month averaged from 2005 to 2008 
at Cariri. In this paper we refer to total wind as the entire 
wind measured by each radar. In each case results represent a 
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Fig. 7. The same as Fig. 5, but averaged from 1998 to 2008 over Cachoeira Paulista. 


monthly mean, with each time/height interval averaged over 
all data available during the month. It can be seen from this 
figure that there is a maximum in the variances around 95 km 
and 02:00 UT, with downward phase propagation in time, 
reaching 85 km around 14:00 UT in February, March, and 
September, indicating a diurnal modulation. It is of interest to 
note that these maxima are centered where the wind changes 
its direction and has the stronger wind shear, represented by 
the dashed line. In the same figure we can see that none of 
the other months exhibit the downward diurnal modulation 
and only show variances increasing with height. 

In Fig. 6, which shows meridional variances at Cariri, we 
can see in the months of February to April, and August to 
October two maxima, the first starting around 89 km close 
to midnight and the second starting around 98/, km at about 
08:00 UT. Both of these show downward phase propagation 
in time and, as in the zonal component, indicate semidiurnal 
modulation. Again, we can see that these modulations ex- 
hibit good correlations with the stronger wind shears. It is 
also interesting to point out here that even from May to July, 


when the winds exhibit strong shears, the variances do not 
maximize in these regions. 

Figure 7 presents the zonal variances (left side) and the to- 
tal wind (right side) over CP averaged from 1998 to 2008. 
Here we note a semi-diurnal modulation with phase descent 
from February to April and in September. The first maxi- 
mum begins at around 95 km at around 02:00 UT and pro- 
gresses down to 80 km at around 12:00 UT, the second be- 
gins at around 10:00 UT at 98 km reaching 80 km at around 
23:00 UT. For these months, the maxima in the variances 
agree with the maxima in the wind shears. On the other 
hand, January and October to December exhibit a diurnal 
modulation with no phase progression and no relation to the 
wind shear. 

Figure 8 shows meridional variances on the left side and 
the total wind on the right, all averaged from 1998 to 2008 
over CP. Here we see strong semidiurnal modulation with 
phase progression from February to April and in Septem- 
ber, again with good agreement between the maxima in the 
variances and wind shears. On the other hand, a diurnal 
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Fig. 8. The same as Fig. 7, but for the meridional component. 

modulation with no phase progression and no relation to 
wind shear is observed from October to December and June 
to July. 

Shown on the left side of Fig. 9 are the zonal variances 
for Santa Maria averaged from 2005 to 2009, with the zonal 
winds on the right. Semidiurnal modulation with phase pro- 
gression is observed from February to April and in October. 
These variations are different from the diurnal variation ob- 
served at Cariri and in agreement with the observations over 
CP. Also, modulations observed in the equinox months are 
similar at the three sites; the maxima of the variances occur 
where the wind shears are largest. Also, diurnal modulation 
is observed in January, November, and December with no 
phase progression and no relation to the wind shear. 

The meridional variances and winds are shown on the left 
and right in Fig. 10, respectively. A semidiurnal modulation 
with phase progression is observed in February and March 
and the maximum of this modulation accompanies the largest 
wind shear. Semidiurnal modulations with phase progression 
are also observed in September and October; however, the 
variance maxima do not follow the largest wind shears. 
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Figure 11 shows the zonal variances and winds for 
SAAMER during 2010. Compared to the observations 
made at Cariri (150 to 400m 2 s~ 2 ), CP, and SM (200 to 
400 m 2 s -2 ) shown in Figs. 5, 7, and 9, respectively, we 
can note that the variances observed at SAAMER (250- 
600 m 2 s -2 ) are larger than those observed at the other three 
locations. Although not so evident as in the zonal compo- 
nent, meridional variances at SAAMER are also larger than 
those observed at the other three sites, comparing Fig. 12 
with Figs. 6, 8, and 10. This is more apparent in a temporal 
average (not shown here), where the SAAMER variances can 
reach values larger than 440 m 2 s -2 while at the other sites 
the maxima were no more than 400 m 2 s -2 . These results are 
consistent with the existence of a hot spot for GW generation 
located in the SAAMER region. The zonal values of the vari- 
ances have similar behavior throughout the year, increasing 
with height. Even though we can observe strong wind shears 
in several months, no associated variance increases are ob- 
served. 
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Fig. 9. Zonal component of variances (left side) and the total wind (right side) averaged from 2005-2009 over Santa Maria. All panels show 
the height/time variation of their respective variable. The altitude is measured in km, the x axis represents Universal Time (UT), and each 
line of the variances contour plot corresponds to 50 m 2 s -2 . The total wind is measured in ms -1 and the dashed line indicates where the 
values are zero. 


Figure 12 shows meridional variances and winds observed 
during 2010 by SAAMER. It appears that a diurnal mod- 
ulation in variance exists from March to May, but with no 
phase progression and no relation to the wind shear. Semid- 
iurnal modulations are observed in January, February, Octo- 
ber (above 88 km), and November, again with no phase pro- 
gression and no relation to wind shear. 

There is also evidence of seasonal variance variations, es- 
pecially at higher southern latitudes. However, the present 
paper focuses only on diurnal and semidiurnal variations. 

4 Discussion 

Diurnal or semidiurnal modulations in GW variances are ob- 
served by all of the meteor radars employed for this study. At 
lower latitudes, Cariri, CP, and SM, these modulations occur 
mainly in equinox months, and with the exception of Cariri 


in the zonal component, all of them are semidiurnal. In addi- 
tion, except for the fall equinox in the meridional component 
at SM, all of the semidiurnal modulations exhibit downward 
time phase progression. Moreover, diurnal modulations oc- 
cur mainly near solstice and, except for the zonal component 
at Cariri, all of the others do not show downward phase pro- 
gression. At SAAMER, these modulations are only observed 
in the meridional component, with diurnal modulation from 
March to May and in the semidiurnal component in January, 
February, October (above 88 km), and November. 

Several studies have shown similar modulation in GW 
variances at other locations. A possible explanation for the 
modulation in the GW fields with tidal periods has been 
given, among other authors, by Isler and Fritts (1996). The 
saturated amplitude for a single GW is typically assumed 
to be |m — c|, where u is the mean wind and c is the phase 
speed of the wave. A medium- or high-frequency GW also 
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Fig. 10. The same as Fig. 9, but for the meridional component. 


is typically refracted by a tidal wind field approximately as 
if it were a local mean wind. For a GW spectrum with phase 
speeds distributed about zero, a diurnal tide in the zonal wind 
can impose either a diurnal or a semidiurnal variation in GW 
variance, depending on the relative isotropy or anisotropy of 
the phase speed distribution. For a diurnal tide, zonal winds 
achieve maxima in a given direction twice a day, suggesting 
a possible semidiurnal variation in GW variance, but with a 
diurnal modulation of propagation direction (e.g., Isler and 
Fritts, 1996). On the other hand, a large diurnal tide superim- 
posed on a strong mean wind can result in a diurnal vari- 
ance modulation due to the asymmetry introduced by the 
large mean wind. Likewise, if the GW spectrum itself is 
anisotropic, the asymmetry of the spectrum will in general 
cause variance enhancements to occur at the same period as 
the lower-frequency motion. These examples, given by Isler 
and Fritts (1996), invoke filtering processes to impart a peri- 
odicity to the GW variance. Good correlation between strong 
increases of GW variances and wind shear, as shown previ- 
ously, are consistent with the idea of tidal filtering of the GW 


spectrum. This is based on the removal of GWs encountering 
tidal shears that decrease their intrinsic phase speeds, caus- 
ing them to approach critical levels and dissipate, and the 
strong amplitude growth of GWs in tidal shear that increase 
their intrinsic phase speeds. Such filtering mechanisms may 
explain some of the modulations at tidal periods and correla- 
tions with tidal wind shears, but they remain to be confirmed 
with more detailed measurements and/or quantitative mod- 
elling accounting for all of the relevant dynamics. 

Another mechanism has been suggested by several au- 
thors, which really represents a sub-set of the dynamics dis- 
cussed above (e.g., Walterscheid, 1981; Thayaparan et al., 
1995; Manson et al., 1998). The mechanism is to again filter 
the GW spectrum, but by critical levels at which a GW has 
a horizontal phase speed equal to the background wind in 
the direction of GW propagation. In reality, however, critical 
levels represent the level below which a GW must dissipate, 
given the saturation limit noted above, according to linear 
theory. 
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Fig. 11. The same as Fig. 9, but over Tierra Del Fuego during 2010. 


A third mechanism associated with diurnal or semidiur- 
nal modulations of potential GW sources, such as deep con- 
vection or mountain wave generation, could enable possible 
explanations for modulations that do not exhibit downward 
phase progression, especially in cases where filtering may be 
a weaker influence due to weak tidal winds or shears. But 
it appears unlikely that these GWs, with either isotropic or 
anisotropic phase speed distributions, would not also exhibit 
downward progression of variance maxima for larger tidal 
winds and wind shears. 

Clemesha and Batista (2008) showed a good correlation 
between wind shear and the wind fluctuation related to GW 
activity. They also suggested that the wind shear in the MLT 
region due to tides could be a source for GWs. It is in- 
teresting to note that as long ago as 1968 Lindzen (1968) 
suggested that breaking of tides could lead to the genera- 
tion instabilities which might produce gravity waves. As is 
well known, the diurnal tide with short vertical wavelength 
is dominant and has a semiannual oscillation with maximum 


amplitude at the equinoxes for the three Brazilian sites (see 
Batista et al., 2004; Andrioli et al., 2009; and Lima et al., 
2007). It is interesting to note that all the variance modula- 
tions that show phase progression and good correlations with 
wind shear occur during equinoxes when the diurnal tide 
has a large amplitude: see, e.g., Fig. 5 (February, March and 
September), Fig. 6 (February to April, August to October), 
Fig. 7 (February to April and September), Fig. 8 (February 
to April and September), Fig. 9 (February to April and Oc- 
tober), and Fig. 10 (February and March). Moreover, GW- 
tidal interactions are strong when tidal amplitudes are large. 
These act both to remove GWs from the spectrum via filter- 
ing, as explained previously, and to excite additional GWs as 
part of the GW dissipation dynamics, though these dynam- 
ics have yet to be fully quantified and understood. As noted 
above, however, there are reasons to expect that GW-tidal 
interactions are a source of additional GWs, as suggested by 
Clemesha and Batista (2008). 
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Fig. 12. The same as Fig. 1 1, but for the meridional component. 


5 Summary and conclusions 

We have presented a modified composite day analysis of all- 
sky meteor radar data extending from 7° S to 53.6° S. We ob- 
served diurnal and semidiurnal modulations of GW variances 
over all the sites. Semidiurnal modulations with downward 
phase progression were observed at lower latitudes, Cariri, 
CP and SM, and they occurred mainly near equinox. On 
the other hand, diurnal modulations occurred mainly near 
solstice and, except for the zonal component at Cariri, did 
not exhibit downward phase progression. At SAAMER at 
53.6° S, these modulations were only observed in the merid- 
ional component, with diurnal modulations from March to 
May and semidiurnal modulations during January, Febru- 
ary, October (above 88 km), and November. After making 
corrections to the variance estimates, as suggested by Andri- 
oli et al. (2013), a number of observed modulations showed 
good correlations with wind shears. The most interesting 


result of this work is that all the cases where the variance 
modulations show phase progression, and good correlations 
with wind shears, occur in the months where the diurnal tide 
achieves large amplitudes, e.g., near equinoxes. 
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